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Abstract 

T\\f6 digital models for estimating daily evaporation and soil moisture 
have been developed as part of the Tellus Project, sponsored by the 
Joint Research Centre of the European Commiinity. One, for grassland 
areas, relates surface temperatures to evaporation, and the other, 
for bare soils, relates day and night temperatures to both evaporation 
and thermal inertia, which may in turn be related to soil moisture 
conteni:. Surface tentperarares may be estinaced using .aeasurements of 
thermal emission, such as those recorded using an infra-red line 
scanner. 

• 

Tliese models have been tested through a flight experiment at Grendon 
Underwood, in Buckinghamshire, bhilst primarily research tools, it 
is hoped t!iat in the future similar models may be used operationally 
to give areal estimates of mean soil moisture and evapotrimspiration. 

To this end, data fresn the Heat Capacity Mapping Mission are being used 
to test the models. w 
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Introduction 

Estimates of soil moisture content and of evapotranspiration are useful 
for hydrological and agricultural putposes; accurate estimates can on 
occasions lead to considerable cost savings » by more effective 
planning of water supply, irrigation and drainage systems. Regional 
estimates of soil moisture content or of evapotranspiration may 
nevertheless be difficult to obtain because of measurement problans. 


•One method of estimating both soil moisture content ai-d svarotrsnsnlrt- 
tion is to use a model of the energy balance of the soil surface and 
of the atmospl^eric boundary layer; such a model may be written in a 
form that most of the input data used are routinely observed, and may 
be adapted so that remotely-sensed measureo'.ents may be included. 


It is possible to estimate the evaporation for one instant by solving 
the transport equation for evaporative heat transfer, and also to 
estimate variables related to instantaneous soil moisture cont^t from 
the diffusion equation. However, the users of soil moisture content 
or evapotranspiration data generally require not merely instantaneous 
values, but cumulative values for particular time intervals, particularly 
daily values. In order to do this, it is necessary either to have a 
time series of measurements of soil moisture or of evapotranspiration, 
or to use a model which estimates these quantities from variables which 
are routinely measured continuously, bhile continuous measurements of 
soil moisture or of evapotranspiration be accurate, it is difficult 
to deploy a large number of such sensors, and so the spatial variations 
of either variable would only bo poorly estimated, although these 
v.-iriaticns may be ccnsideriblj. Ti.e e:''plo',.T.';nt of a digital r.JO'del is 
thus very dcs irabl *. In :rdvM* cr.-jl:/ nodoLs, it is noce>iar/ 
to estimate t’ne components oc the cncr?>' balance of the soil surface 


for each of the ttoe steps of measurements input. 

The txvorzY balance of the soil or crop surface is governed by a 
emtinuity equation 

G + H <*• L.E - 0 

Vrliere is the net radiation flux, G the soil heat flux, the 
sensible heat flux in the atmospheric boundary layer and L.E. the latent 
heat flux. L.E is tlie product of E, the evapotranspiration flux, 
and L, the latent heat of vaporisation of water per unit mass. 

Caiventionally, dcnviiward fluxes are given a positive value. 

COTsider the componaits of the radiation balance in more detail. The 
net radiation flux is the sura of the inconing and outgoing short and 
long wave radiation fluxes: 

- (1 . ag)Rg + (1 - a^)K^ - c^oT^ (2) 

^.ere is tlie measuring sliort wave and the incoming long wave 
radiation, and o, die short and long-wave soil/crop reflection 
coefficients, e^the soil/crop emission coefficient, o the Stefan 
Botrmium constant, and T, the soil/crop surface temperature. Radiation is 
considered short-wave up to a wavelength of l.l ym. For long-wave radiation, 
e^* 1 ■ so equation (2) may be ivritten 

R,, • (1 - a)s Rj + 

The soil heat flux G throu^jh any honeon ; ■ is related to the 
temperature gradient in the soil 5T/3i and the thermal conductivity X, 
both evaluatoJ at* 
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This relationship leads to a form of the diffusion equation ' '“ 

" \ • 

■* ' I = 

'* 3 r^ . r ^ 

5; (xjj) ■ <v 55 

where C^ is the volumetric heat capacity of the soil. This equation may be 
solved numerically for each level z' to gi\-e a temperature pi*ofile using a 

finite eicanenr. or finite tiiffe fence i.iethod, takLig iha :»oil/crop jurf'-uai 
temperature as the upper boundary condition, and sane known temperature 
at a given depth in the soil as a lotpfer boundary condition. 


The sensible heat flux may be estimated using a transport equation, 
depending on the air temperature T^, the soil/crop surface, T^., the 
aerodynamic resistance, r , and the air density p and specific heat c 


H ■ pc. 




The latent heat flux Jiuiy also be expressed as a transport equation 


K - 

h * P's 


where y is the psychroraetric constant, e. the water vapour pressures 

£l 

in the atmosphere layer, e the saturated water vapour pressure at 
temperature T^, and r^ the stomatal diffusion resistance to water 
vapour transport. Tiie quantities and may be estimated from 
and T . 


13)' coi-^parin? rt), (-H , (5) and f61 , it may be seen that the components 
of r!:-> -l- ,;iav be o.^nresseJ as fanaticas o; the air 


temperature, tiie temperature at a given depth in the soil, and the 

soil surface temperature, as well as ga the incoming radiation and 

^ . the aerodynamic and stomatal resistances. Most of these variables may 

be measured or estimated from routine observations, which are widely 

collected continuously in time. Tlie exception is the crop/soil surface 

* 

• temperature, T , which is not routinely measured except at a few sites. 

ttewever,. if a time series of estimates of the crop surface tenperaturcs 
K may be made then the evapotranspiration flux E and the soil moisture 

may be estiniHCcd. Such estimates may be provided using digital models. 

, f 

Digital models 

Equations (3), (4), (5), and (6) may be re-written to give an expression 
. for the soil/crop surface temperature in teims of routine mwsurements. 

The initial estimate of the surface temperature will probably not 
satisfy the continuity requirement of equation (1), because of 
sijiplifications which arc necessary to estimate parameters such as r , 
and so it is necessary to use an iterative numerical optimization 
technique to ensure that the condition in equation (1) is satisfied. 

Tl\is procedure may be mlojited to give estimates of the surface 
temperature corresponding to every tune step of a set of routine 
measurements of air temperature, wet-bulb temperature, wind speed and 
incoming short and Icng wav? radiation. 



Althoug!: continuous measurements of surface temperature at not made 
routinely, it is nevertheless possible to make instantanteous 
luoasurcmonts, by using a radiometer to measure the theiral emission 


from tliO soil surface. Tr.c measured emission may be related to the 
suvtnco tcr-'orafiro, ns snov.n Ln the ne.xt sccti 


vH. It is possible to 
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compaxe.the measured surface temperatures with those estimated from a 

B 

model, and then to adjust the parameters of the model until the 

la 

estimated surface temperatures are the sane as the measured sur&ce 
teiq>eraturos at the same time. These parameters such as the aero- 
• ^mamic resistance r^^ and the thermal heat capacity of the soil Cy, • 

may then be used to give more accurate estimates of daily evapo- 
■ transpiration or soil moisture. Two models luve been developed 

using these methods. Gne model, the "Tergra" model, developed by 
Soer (1^77), is for use in grasslnixd areas vuid, using one ’.noasuremenc 
of the soil surface tenperature, gives estimatos of the daily evaporation. 
The other, the 'Tellus" model, developed by Rosena (Rosema at al., 1978), 
is for use on bare soil, and estimates both daily evaporation and 
thermal inertia usuig measurements of botii day-time and night-time 
surface temperatures. The thermal inertia p may be defined as 
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P - (XCy)i 

The relationship between thermal inertia and soil moisture may be 
determined cAperimentally. Both models are layered miodels, which 
require known boundary conditions at a given height in the atmospheric 
boundary layer .and at a given depth in the soil in order to solve the 
equation (1) using equations (3) , (4) , (5) and (6) . A time series 
of measurements of incoming radiation, temperature, vapour pressure, 
and wind speed at (say) 2 m is then input to both models, together 
with soil conditions at a sufficient depth in the soil that they may be 
assumed not to vary during the period over which the models are applied. 
For each forw;ird tine step, equation (1) is solved to estimate the 
coiiipcr.cnts of t;;o energy balance, which may then be used to estimate 


the quantities of 


interest Ln each model. 


Uacii ::;cdel will now he 



considered briefly, to point to the differences in the way equation (1)' 
is solved in each case. 

The Tergra model for grassland areas uses a Businger-Dyer approach to 
estimate the turbulent diffusion resistance r^, which is used in 
equations (S) and (6) . This gives r^ as a function of wind velocity, 
the stal)llity of the atmosi^eric boundary layer just above the surface, 
and the nature of the surface. Different conditions are applied under 
stable, noutnl and unstable act^ospheric conditions. Hie stcmatal 
resistance, r , used in equation (6), is found from a set of 
ei!q>irical relationships wlildi relate -it to the crop height, the 
short-wave radiation and the leaf-wat^'r pressure (Soer, 1977) . 

The soil heat flux is estimated using an explicit finite difference 
method with the spacings between each node being equal. Tlie iteration 
of the surface touperaturc is also simple, the surface temperature being 
altered by successive small iimounts until an energy balance has been 
achi«\red. Tills approach is satisfactory where the components of the 
energy balance are not changing rapidly between successive time steps, 
so that the uiitial estimate of the surface temperature, (the teiiiperature 
estimated firm tlic previous time step) is close to the final estimate. 

As noted ab«5v'e, the Tergra model directly estimates the daily 
evaporation. Oianges in the soil moisture may also be estimated from 
the evapotranspiration fliuc, E, estimated in equation (6). 

Given the cvapornti.cn, soil moisture is chtalr.ed by a set of rclaticnsh.tps 
as follo’.-.s. For "rnssland, E na'* bo e-mressed ns 






where is the leaf water pressure, the soil water pressure, 
^plant resistance for water transport, and the soil 

hydraulic resistance. Given estimates of and the resistances 
*p?ant *^soil* ^ estimated. These may be found using 

enq)irical relationships with other variables (Soer, 1977) and then 
may be related to the volumetric water content, 6, by 



S-S, 

rrs: 


( 8 ) 


idiere S is the saturation, defined as 6/0^ (6^ being the moisture 
content at saturation) . is the rest saturation, m is a pore size 

distribution factor, and is a re-scaled soil water pressure 

equal to where is the air entry value. 

The Tellus model for bare soil sites uses titfo separate estimates of 
surface temperature, ideally at the times of maximum and minimum, 
themal emission (about 02C0 and 1400 local time) , and so allows 
two unknowns to be estimated directly, without the need to use 
empirical relationships such os (7) and (S) to derive one of these 
unknowns. In this model, themal inertia and evaporation and 
trunspiratica are estimated; an experimentally-derived relationship 
may be found between themal inertia and soil moisture, hhilst 
the general structure of the model is sLailar to that used in • 
the Tergra model, there are some difflftrences. The soil heat flux is 
again found using an explicit finite difference method, but in the 
Tellus model the spacing bet«.veen nodes Is allowed to increase dotvnv;ards, 
using a Dufort .uid Frankel method. The iteration at each time step 
is also rather faster than for the Tergra model, as a Nev/ten -Ralph son 
methed is this converges rapidly. 


<» 


H» most criticol input . , in ths sense of being roost difficult 

to estiroete, is the surface temperature. Other inputs are either 
measured, or may be estimated. The roughness length, which may be 
difficult to estimate end is used in the estimation of the aerocfynamic 
resistance, is estimated using empirical relationship:), as suggested 
by Monteith (1973) . 


Esthnation of surface temperature using radionetew 

Radiometers me.nsure the radiant energy received at the sensor in a 
set of wavelengths. IVhen directed at an area of the ground surface, 
the energy received must first be related to tlie radiant emission 
from tliat area which in turn roust be related to the teo^erature of 
that area. The radiometer itself is calibrated using two black ' 
body radiators, one hotter and one colder than the area of ground 

t 

being studied. 

The relatiori^lp between emission and tenperature usually employed - 

Stefan's I.aw - is only valid when the emission is integrated over 

all wavelengths. However, it has been shown (eg Scarpace, 1974) 

that to use Stcion's ktw for the 3 - 14 um tvaveband often adopted 

for radiometers introduces only very small errors, usually less 

than 0.03K. The temperature of the soil surface nay then be estimated^ 
from the measurements E* by a relationship of the form 

E. - (1 - e^) CjOTj ♦ e^cTj (9) 


where is the cmissivity of the crop or soil surface, is^ the 
sky cmissivity, T is the apparent sky tcr.:pefature (K) T is the crop 

9 W 


surface temperature, uiiJ c is Stcfan-Bolczirunn's constant. If tiio 


emission as measured at the railicT.eter 
into e<ir.al-sLroJ ^teps between the r..o 


is digitised by dividing 
C’.libratioa black-bcay 


it 


omissions on the radiometer, the measured emission E» may be estimated; 
if (for example) the data are digitized into 8>bit nunbers; scanner step 
2SS vdll be set equal to the emission at the upper calibi-ation 
ten^vature, scanner step 0 set equal to that at the la- or calibration 
temperature, and E« is estimated as 

255 X c,eT*’ ♦ c*cT^g- x SS - c*oTX x SS 

E. - 2 0 (jOj 

where e« is t!ic cmissivity of the "black body" calibration on the 
I’adi-iinecw^r, ;a*.d ^255*’*’’* twniperaturis of the Icwi.-r and upper calibration 
"black bodies", and SS tlte scanner stop equivaleiit to the measured ground 
temperature. By substitution of (10) into (9), the soil surface temperature, 
T^, may be found izi terns of the scanner step, It is necessary to 
know the eraissivlty of the soil surface, and the apparent sky temperature 
and its emissii^ity, as wll as the emissivJt>’ of the calibration black 
bodies, tlieir temperatures and the scanner step of the area being examined. 
These may be found froci s;uii()]e ground measurements: the tenq)erature 
estimated is not very sensitive to those parameters for typical bare soil 
or short vegetation, where tlic cmissivity is high. Figures 1 rid 2 show 
tlie absolute errors which aro^’pt jduced for a typical set of values, if •' 
the sky tcmpenitnro and crop tinissivity are incorrectly estimated. For 
this example, the correct values are a sky temperature of 2"^ K, croo and sky 
emissivitiesof 0.95, a black b^xly emissivity of 0.99 ttid calibration 
temperatures of 273 K and 291 K and a crop temperacore of 281 K. The 
sky tenperature was alterei by * 20K in IK steps and the crop emissivity 
was allowed to vaiy from 0.76 to 1.0, the ihsolute errors in 
temperature being noted. It nay he seen chat the sky temperature need only be 
found withi:\ 13K of th.e correct value, and the emissivity need only be 
estimated wiclun 0.C5 for the- cr'?p .'Ui^acc te;r,peratut'e to be estircitcd to 


within O.SK of tht correct value. 


Atmospheric effects are also very in^rtant, particularly vhen 
measurements from satellites are used. Although It is possible to 
construct eqiuitions to estimate the absoxption, scattering and 
re-emission of radiation using Mie theory, it is not practically 
possible to collect enough data to allow these equations to be solved, 
and so approximations are necessary. One such approximation where the 
a5si<a|*tiui\ Is rode tiut ail trouiacian is absorbed und ro-e>aitted by 
water vapour has been suggested byRangaswaK^sr al. (1978). In all 
cases it is necossaiy to have information about ground tenq>eratures and 
about the temperature and humidity profile of the atmosphere. 

Application of the models 

This project is part of the Tollus project, sponsored by the Joint 
Research Centre, Ispra. As part of this project a Joint Flight 
Experiment was perfonned at Crendon Undervood, BucUinghamshire in 
Septenfeer 1977. Most of tht^, couq^onents of the energy balance were 
measured dircctiyi^ap.irt from actual evaporation and transpiration, 
which are difficult to mcaiure. Two sites were monitored in detail, 
one field of hare soil and one field of grass of length S - 10 an. 

Two flights were made, using a Daedelus DS-12S0 multi-spectral 
Scanner with a thermal infra-red sensor in the S - U un wavebend, 
one flight being at 5 am and the other at 2 pm BST. Simultaneous greund 
neasi-.rcmcnts with PKT-5 radiometers ivere made; these treasure emission 
in the S - 14 um *.<aveband with a 2^ field of view. Temperatures 
estitutod from the aircraft and ground-mounted radiometers were ver>' 
.sunilar for bot!i day-ture and night-time iat;i. 


prelimijuiry I'esults of moUelling are presented in Table 1. Nfeon 

soil surface temperatures were estimated fron radiometer data for the bare 

soil site and the graasland site. Both morning and afternoon temperature 

estimates were input to the Tellus model » and the morning temperature 

estimate only was input to the Tergra model. For comparison^ the mean. 

wlumetric soil moistures are included, ta^ from 15 on core sanples 

for each site. ‘Hic actual evaporation for 13tli September 1977 is difficult 

to estimate, flowever, if the lun-off from the River Ray catcliment 

\nthin ’.vhich the test areas are sited is si^tricted frtn the rainfall 

then some ide.a of the evaporation may be obtained if change in water 

storage \d.thin the catchment is neglected. This obviously involves 

many assuiptions about the behaviour of the catchment, but if this 

exercise is performed for data for the month of SeptCTsber 1977, and then 
# 

divided to give a daily figure, then this may give an estimate of the 
.actual evaporation to an order of magnitude. Such a figure is included 
in Table 1. As 13th Septeuiier, the day of the flight experiment, was 
clear and rainless, the actual evaporation and transpiration were probably 
above the mean figure given. Ihe Tergra model rather underestimated the 
evaporation and transpiration, but as the results are only preliminary 
the reasons for this are not yet understood. 

Further analysis of these data is continuing. It is hoped that data 
from the Heat Capacity Mapping Mission for the U.K. area will also be 
used as inputs to tne models. Much further work is needed to test the 
models outlined here. Nevertheless, the preliminar/ results presented 
are very encouraging and so it may be hoped that regional estimates 
of soil moisture .aiid of evaporation and transpiration will soon be 
available using this approach. 
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